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Abstract — The design of modulation schemes for the physical 
layer network-coded two way relaying scenario was studied in 
|1 1, 1 3 1, |4| and |5|. In |7| it was shown that every network coding 
map that satisfies the exclusive law is representable by a Latin 
Square and conversely, and this relationship can be used to get 
the network coding maps satisfying the exclusive law. But, only 
the scenario in which the end nodes use M-PSK signal sets is 
addressed in |7| and |8|. In this paper, we address the case in 
which the end nodes use A/-QAM signal sets. In a fading scenario, 
for certain channel conditions ye* , termed singular fade states, 
the MA phase performance is greatly reduced. By formulating a 
procedure for finding the exact number of singular fade states for 
QAM, we show that square QAM signal sets give lesser number 
of singular fade states compared to PSK signal sets. This results 
in superior performance of A/-QAM over M-PSK. It is shown 
that the criterion for partitioning the complex plane, for the 
purpose of using a particular network code for a particular fade 
state, is different from that used for M-PSK. Using a modified 
criterion, we describe a procedure to analytically partition the 
complex plane representing the channel condition. We show that 
when Af-QAM (M > 4) signal set is used, the conventional XOR 
network mapping fails to remove the ill effects of ye je — 1, which 
is a singular fade state for all signal sets of arbitrary size. We 
show that a doubly block circulant Latin Square removes this 
singular fade state for M-QAM. 

I. Preliminaries and Background 

We consider the two-way wireless relaying scenario shown 
in Fig.l, where bi-directional data transfer takes place between 
the nodes A and B with the help of the relay R. It is 
assumed that all the three nodes operate in half-duplex mode. 
The relaying protocol consists of the following two phases: 
the multiple access (MA) phase, during which A and B 
simultaneously transmit to R using identical square M-QAM 
signal sets and the broadcast (BC) phase during which R 
transmits to A and B using either a square M-QAM signal 
set or a constellation of size more than M. Network coding 
is employed at R in such a way that A (B) can decode the 
message of B (A), given that A (B) knows its own message. 

A. Background 

The concept of physical layer network coding has attracted a 
lot of attention in recent times. The idea of physical layer net- 
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(b) BC Phase 
Fig. 1. The Two Way Relay Channel 

work coding for the two way relay channel was first introduced 
in Ql, where the multiple access interference occurring at the 
relay was exploited so that the communication between the 
end nodes can be done using a two stage protocol. Information 
theoretic studies for the physical layer network coding scenario 
were reported in j2), J3}. The design principles governing the 
choice of modulation schemes to be used at the nodes for 
uncoded transmission were studied in |4). An extension for 
the case when the nodes use convolutional codes was done 
in (5j. A multi -level coding scheme for the two-way relaying 
scenario was proposed in [6|. 

It was observed in [4| that for uncoded transmission, the 
network coding map used at the relay needs to be changed 
adaptively according to the channel fade coefficients, in order 
to minimize the impact of the multiple access interference. 
The Latin Square scheme was studied in |7J, |8) for two way 
relaying using Af-PSK signal sets at the end nodes. 

B. Signal Model 

Multiple Access (MA) Phase: Let S denote the square M- 
QAM constellation used at A and B, where M = 2 2A , A 
being a positive integer. Assume that A (B) wants to transmit 
a 2A-bit binary tuple to B (A). Let v : F 2 2a — > S denote the 
mapping from bits to complex symbols used at A and B. Let 
v{sa) = xa, v(sb) — xb G S denote the complex symbols 
transmitted by A and B respectively, where sa,sb € F 2 2a. 



The received signal at R is given by, 

Vr = h A x A + h B x B + z R , 

where h A and h B are the fading coefficients associated with 
the A-R and B-R links respectively. The additive noise zr 
is assumed to be CAf(0, u 2 ), which denotes the circularly 
symmetric complex Gaussian random variable with variance 
a 2 . We assume a block fading scenario, with the ratio h B /h,A 
denoted as z = jeJ , where 7 € K + and — tt < 9 < n, is 
referred as the fade state and for simplicity, also denoted by 

Let 5r(7, 8) denote the effective constellation at the relay 
during the MA Phase and d m i n {l& ) denote the minimum 
distance between the points in Sn(-f, 9), i.e., 

= {xi + je 3e x j \x l ,x j e5}, 
dminhe 30 ) = min ; | (x A - x' A ) + je 30 (x B - x' B ) \. 

(x A ,x B ),(x' A ,x' B )£S 
{x A ,x B )=jt(x' A ,x' B ) 

(1) 

From ([TJ, it is clear that there exists values of je J for 
which d min (je 30 ) = 0. Let H = {je 30 G C\d min (j, 9) = 0}. 
The elements of T~L are said to be the singular fade states (7). 
The set W depends on the signal set used. For example when 
je 3 = (l+j)/2, the effective constellation Sr(j, 9) has only 
12 (<16) points when 4-QAM signal set is used at nodes A 
and B. Hence (1 + j)/2 G U for 4-QAM. 

Let (£ A ,x B ) G S 2 denote the Maximum Likelihood (ML) 
estimate of (x A ,x B ) at R based on the received complex 
number ur, i.e., 

{x Al x B ) = argmin \y R - h A x' A - h B x' B \. (2) 

(x y A y B )es 2 

Broadcast (BC) Phase: Depending on the value of -fe 3 , R 
chooses a map M 1,0 : S 2 — > S\ where S' is the signal set (of 
size between M and M 2 ) used by R during BC phase. The 
elements in S 2 which are mapped on to the same complex 
number in S' by the map J\4 7,0 are said to form a cluster. 
Let {£ 1; £ 2 , •••) £1} denote the set of all such clusters. The 
formation of clusters is called clustering, and the set of all 
clusters is denoted by C 7,0 to indicate that it is a function of 
je . The received signals at A and B during the BC phase 
are respectively given by, 

Va = h' A x R + z A , y B = h' B x R + z B , (3) 

where xr = M. 1 ' (£ A ,x B ) G S' is the complex number 
transmitted by R. The fading coefficients corresponding to the 
R-A and R-B links are denoted by h' A and h' B respectively 
and the additive noises z A and z B are CAf(0, a 2 ). 

In order to ensure that A (B) is able to decode B's (A's) 
message, the clustering C' e should satisfy the exclusive law 
@, i-e., 

M~<' e (x A ,x B ) ^ M~<< (x' A ,x B ), lorx A ^x> A ,Vx B ES, 1 
M~<' (x A ,x B ) ^ M~ f ' e (x A ,x' B ), forx s ^ x' B ,\/x A 6 5. J 

(4) 



The cluster distance between a pair of clusters £j and Cj 
is the minimum among all the distances calculated between 
the points x A + ^e 36 x Bl x' A + ^ye^x'g G Sr(j,9) where 
(x A ,x B ) G Ci and (x' A ,x B ) G Cj |7j. The minimum cluster 
distance of the clustering C is the minimum among all the 
cluster distances, i.e., 

<Cm(7e J ' e )= min | (x A - x' A ) +ie je (x B - x' B ) \. 
es 2 , 

M~i- e { XA ,x B )^M~<' e (x' A ,x' B ) 

The minimum cluster distance determines the performance 
during the MA phase of relaying. The performance during the 
BC phase is determined by the minimum distance of the signal 
set S' . For values of je 1 in the neighbourhood of the singular 
fade states, the value of dmin^e 30 ) is greatly reduced, a 
phenomenon referred as distance shortening. To avoid distance 
shortening, for each singular fade state, a clustering needs to be 
chosen such that the minimum cluster distance at the singular 
fade state is non-zero and is also maximized. 

A clustering C is said to remove a singular fade state h G 
T-L, if d%ni n {h) > 0. For a singular fade state h G W, let 
denote a clustering which removes the singular fade state 
h (if there are multiple clusterings which remove the same 
singular fade state h, consider a clustering which maximizes 
the minimum cluster distance). Let Cr = {C^ h ' : h G H] 
denote the set of all such clusterings. Let d m in(C^ h \ 7', 9') 
be defined as, 

d mm (C {h \l',e')= min \(x A -x' A )+y'e i$ ' (x B -x' B )\. 

(x A ,x B ),(x A ,x B ) 

es 2 , 

M< h Hx A ,x B )^M< h hx' A ,x' B ) 

The quantity d m i n (C^ h \j,' 9') is referred to as the mini- 
mum cluster distance of the clustering evaluated at Ye 3 ® . 

In practice, the channel fade state need not be a singular 
fade state. In such a scenario, among all the clusterings which 
remove the singular fade states, the one which maximizes 
the minimum cluster distance is chosen. In other words, for 
^'e 30 ^ W, the clustering C 7 ,e is chosen to be C* h >, which 
satisfies d mm {C {h \i ',#') > <WC {h '>,V, 6'), V/i ^ h' E 
H. Since the clusterings which remove the singular fade states 
are known to all the three nodes and are finite in number, the 
clustering used for a particular realization of the fade state can 
be indicated by R to A and B using overhead bits. 

Example 1: In the case of BPSK, if channel condition is 
7 = 1 and 9 = the distance between the pairs (0, 1)(1,0) 
is zero as in FigQa). The following clustering removes this 
singular fade state. 

{{(0,1)(1,0)},{(1,1)(0,0)}} 

The minimum cluster distance is non zero for this clustering. 
The contributions and organization of the paper are as follows: 

• A procedure to obtain the number of singular fade states 
for PAM and QAM signal sets is presented. 

• It is shown that for the same number of signal points M , 
the number of singular fade states for square Af-QAM 



(fO Effective Constellation at the relay for 7=1 and 9 = 



(It) Effective Constellation at the relay for 7 = 1 and 9 = 7 



Fig. 2. Effective Constellation at the relay for singular fade states, when the 
end nodes use BPSK constellation. 



is lesser than the number of singular fade states for M- 
PSK. The advantages of this are two fold - QAM offers 
better distance performance in the MA Phase and QAM 
requires lesser number of Latin squares (i.e., a reduction 
in number of overhead bits). 

• It is known from fl"5) that the removal of the singular fade 
state z=l assumes greater significance in a Rician fading 
scenario. The bit-wise XOR map removes this singular 
fade state when A/-PSK signal set is used at nodes A 
and B for any M. It is shown that XOR mapping cannot 
remove this singular fade state for any A/ -QAM and a 
different mapping is obtained to remove this singular fade 
state for M -QAM. 

• Inspired from [8|, the problem of partitioning the entire 
complex plane into clustering independent region as well 
as clustering dependent region is addressed. The approach 
followed for Af-QAM signal set needs to be different 
from that used for M-PSK in (8). 

• The region associated with each singular fade state in 
the complex plane is obtained analytically for M-QAM 
signal set used at nodes A and B. This helps in associating 
a Latin Square corresponding to that singular fade state to 
the said region like in A/-PSK. By simulation it is shown 
that the choice of 16-QAM leads to better performance 
for both the Rayleigh and the Rician fading scenario, 
compared to 16-PSK. 

The remaining content is organized as follows: 

In Section [II] we discuss the relationship between singular 
fade states and difference constellation of the signal sets used 
by the end nodes. We present expressions to get the number of 
singular fade states for PAM and square QAM signal sets in 

it 



Subsections II-A and II-B respectively. In Subsection II-C 



is proved that the number of singular fade states for M-QAM 
is always lesser than that of Af-PSK signal sets. In Section 
III the clustering for singular fade states is obtained through 
completing Latin Squares and a Latin Square for removing 
the singular fade state z = 1 is analytically obtained for PAM 



and QAM signal sets. In Section IV channel quantization for 
M-QAM signal set is discussed. In particular, the channel 
quantization for the entire complex plane when nodes A and B 
use 16-QAM signal set is obtained. In Section [V] simulation 
results are provided to show the advantage of Latin Square 
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(a) y/M PAM constellation 



A5 = -2(v / M-l) + 2n 



Fig. 3. 



(b) Difference Constellation 



M PAM constellation and difference constellation for v M 



scheme for QAM over XOR network coding scheme as well 
as Latin Square scheme for PSK signal sets under Rayleigh 
and Rician fading channel assumptions. 

II. Singular Fade states and Difference 
Constellations 

The location of singular fade states in the complex plane 
for any constellation used at end nodes can be characterised 
in the following way. If node A uses a constellation S\ of size 
Mi and node B, a constellation S2 of size M2, the singular 
fade states z = je^ e are of the form 



z = je J 



(5) 



x B - X B 

and is obtained by equating xa + je^xs and x' A + ^e^ B x' B 
for xa,x a G Si and xb,x' b G £2- Henceforth, throughout 
the paper, we assume both the end nodes use the same 
constellation, S. 

A. Singular Fade States of PAM signal sets 

Consider the symmetric %/M-PAM signal set given by 

S = {-(\/M-l) + 2n}, n6(0,-,v / M-l). 

The difference constellation of S is 

AS = {x-x' : x, x' G S} 

and can be written in the form 

AS = |-2(%/M - 1) + 2n| , n G (0,- •• , 2(VM - 1)). 

For example, the 4-PAM signal set and it's difference constel- 
lation are given in Fig |3(a)| and Fig |3(b)| respectively. For each 
of the difference constellation points, the pair in the signal 
set which corresponds to this point is also shown. We will 
often consider only the first quadrant of AS, denoted as AS + , 
which for a general complex signal set is given by 

AS + = {a : real(a) > 0, imaginary(a) > 0}. 

The following lemma gives the number of singular fade states 
for PAM signal sets. 

Lemma 1: The number of singular fade states, for a regular 
a/M-PAM signal set, denoted by N^^j_ pAM ^ is given by 



(Vm-pam) 



VM-1 

4 E n 

n=l 



n 

p\n 



1 



(6) 



where p\n stands for prime number p dividing n. 

Proof: See Appendix |A| ■ 
Example 2: Consider the case of 4-PAM (M = 16) signal 
set as given in Fig|3] There are 2(\/M— 1) = 6 non-zero signal 
points in the difference constellation. Scaled AS + has (vm- 
1) = 3 signal points-{l, 2, 3}. And there are 14 singular fade 
states- 

1 1 2 9 ■? 3 _1 _1 ~ 2 9 ~ 3 
' 2' 3' 3' ' ' 2'~ '~2~'T'T' _ ' '~2~ 

Calculating |6]) also gives -/V^p^,/) = 14- Calculating 
similarly, we find that for 8-PAM signal set, there exists 70 
singular fade states. 

B. Singular Fade States for QAM signal sets 

We consider square M-QAM signal set S = {A m i + 
jA m Q}, where A m j and A m Q take values from the \/M-PAM 
signal set -(y/M - 1) + 2n, n £ (0, • • • , y/M - 1). We 
use the bijective mapping p : S — > Zm = {0, 1, • • • , M — 1} 
given by 

A mI +jA mQ -> i[(%/M-l+A m/ )%/M+(VM-l+^ mQ )] (7) 

for concreteness, even though our analysis and results hold 
for any map. The difference constellation AS of square M- 
QAM signal set forms a part of scaled integer lattice with 
(2y/M - l) 2 points. The 4-QAM signal set with the above 



. 3 



mapping and its difference constellation are shown in Fig 4(a) 



and in Fig 4(b) 



In a practical scenario, there can be an average energy 
constraint E to be satisfied at nodes A and B. In such a 
case, we use a scaled version of the M-QAM signal set 
given by -^{A m i + jA m Q}, where p is chosen so as to 
meet the constraint E. As a special case, for E = 1 (unit 
normalisation), p = gpfzrn ■ ^ ma y be notec l nere that the 
values of the singular fade states (for a particular choice of 
S) are unaffected by the energy constraint E. However, the 
minimum distances of the constellation S and the effective 
constellation <Sr(7, &) are dependent on the choice of E. In 
particular, for unit normalisation, 



d m in (M-QAM) = 



6 



(8) 



VP V M - 1 
Compare this to the case with A/-PSK whose 
d mm = 2sin(7r/M). For M = 16, d mm (16-QAM) = ^| 
> ^mm(16-PSK) = 2sin(7r/16). This has a detrimental effect 
in the performance during the MAC phase. 

The signal points in the difference constellation of M-QAM 
are Gaussian integers scaled by J~p. To get the number of 
singular fade states for square QAM signal sets, the notion of 
primes and relatively primes in the set of Gaussian integers is 
useful. 

Definition 1: fTJ) The Gaussian integers are the elements 
of the set Z[j] = {a + bj : a,b € Z}, where Z denotes the set 
of integers. A Gaussian integer a is called a Gaussian prime if 
the Gaussian integers that divide a are: —j, a, —a, aj 

and —aj. The Gaussian integers which are invertible in Z\j] 



0. ! . 2 

(a) 4— QAM constellation 

! A5+. 



(b) The Difference Constellation 



-1.5 -1 -0.5 0.5 1 1.5 

(c) Singular fade states 

Fig. 4. 4— QAM constellation, its difference constellation and singular fade 
states 



are called units in Z[j] and they are ±1 and ±j. Let a,j3 £ 
If the only common divisors of a and j3 are units, we 
say a and j3 are relatively primes. 

Lemma 2: The number of singular fade states for the square 
M-QAM signal set, denoted by Nm-qam is given by 



N M - 



M—QAM 



80(A5+) 



where cj)(AS + ) is the number of relative prime pairs in AS + . 

Proof: All the possible ratios of elements from AS give 
singular fade states. We consider only ratios in AS + and 
multiply the number of such possible ratios by a factor of 4 
to account the ratios with points in all the other quadrants. To 
avoid multiplicity while counting, we take only relative prime 
pairs in AS + and every such pair (a, b) gives two singular fade 
states a/b and b/a. Because of this, the multiplication factor 
becomes 8. Finally, the constant term 4 is added to count the 



units. 



Example 3: For 4-QAM signal set shown in Fig |4(a)| the 
number of singular fade states, N^-qam is 12. Scaled AS + 
has only two elements {1, 1 + j} in this case, as shown in 



Fig 4(b) They form a relatively prime pair. The singular fade 



4(c) 



states ±1, ±j, ±1 ± j, ±l±j are shown in Fig. 

Example 4: Consider the case of 16-QAM signal set. It can 
be verified that there are 48 distinct pairs of relative primes, 
and from Lemma [2] Nie-QAM turns out to be 388. 

TABLE I 

Comparison between M-PSK and M-QAM on number of 

SINGULAR FADE STATES 



M 


A/-PSK 


M -QAM 


4 


12 


12 


16 


912 


388 


64 


63,552 


8388 



C. Singular fade states of M-PSK and M-QAM signal sets 

In this section we show that the number of singular fade 
states for A/-QAM signal sets is lesser than that of M- 
PSK signal sets. The advantages of this are two fold- QAM 
offers better distance performance during the MA phase and it 
requires lesser number of overhead bits during the BC phase, 
since the required number of relay clusterings is lesser in the 
case of QAM compared to PSK. 

Lemma 3: The number of singular fade states for M-QAM 
signal set is upper bounded by 4(?i 2 — n + 1), where n = 
i[(2\/M- l) 2 - 1], which is same as 4M 2 - (2M- 1)VM+ 
!)• 

Proof: There are {(2\/M~— l) 2 — 1] non zero signal points 
in AS* which are distributed equally in each quadrant, i.e., 
the number of signal points in AS + is |[(2vM- l) 2 — 1] 
which we denote by n. The maximum number of relatively 
prime pairs in a set of n Gaussian integers is "( n ~ 1 ) , Since 



an upper bound is of interest we substitute this in Lemma [2] 
instead of <j>(AS + ). This completes the proof. ■ 

In JS), it is shown that the number of singular fade states 
for M-PSK signal set is M(^ - 4f + 1), of 0(M 3 ). From 
Lemma [3] an upper bound on the number of singular fade 
states for A/-QAM is of C(M 2 ). Hence, the number of 
singular fade states for M-QAM signal set is lesser than that 
of M-PSK signal sets. 

The advantage of using square QAM constellation is more 
significant in higher order constellations as shown in Table [I] 
For instance, 64-QAM has 8,388 singular fade states where as 
64-PSK has 63,552 singular fade states and the relay has to 
adaptively use 63,552 clusterings. So, with the use of square 
QAM constellations the complexity is enormously reduced. 

III. Exclusive Law and Latin Squares 

Definition 2: flO) A Latin Square L of order M with the 
symbols from the set Z t = {0, 1, • • • , t — 1} is an M x M 
array, in which each cell contains one symbol and each symbol 
occurs at most once in each row and column. 



In J7], it is shown that when the end nodes use signal sets of 
the same size, all the relay clusterings which satisfy exclusive- 
law are equivalently representable by Latin Squares, with the 
rows (columns) indexed by the constellation point indices used 
by node A (B) and the clusterings are obtained by placing into 
the same cluster all the row-column pairs which are mapped 
to the same symbol in the Latin Square. 



A. Removing Singular fade states and Constrained Latin 
Squares 

The minimum size of the constellations needed in the BC 
phase is M, but it is observed that in some cases relay may 
not be able to remove the singular fade states with t = M 
and t > M results in severe performance degradation in the 
MA phase p). Let (k,l) and (k',l') be the pairs which give 
the same point in the effective constellation Sr at the relay 
for a singular fade state, where k, k', I, I' E {0, 1, M — 1} 
and k, kl are the constellation points used by node A and 
1,1' are the corresponding constellation points used by node 

B. If these are not clustered together, the minimum cluster 
distance will be zero. To avoid this, such pairs should be in 
the same cluster. This requirement is termed as singularity- 
removal constraint |7 j. So, we need to obtain Latin Squares 
which can remove singular fade states and with minimum 
value for t. Towards this end, for a given singular fade state 
z E W, initially we fill the M x M array in such a way that 
the slots corresponding to a singularity-removal constraint are 
filled using the same element. Similarly, we fill in elements for 
the other singularity removal constraints for the given singular 
fade state. This removes that particular singular fade state. 
Such a partially filled Latin Square is called a Constrained 
Partially-filled Latin Square (CPLS). After this, to make this 
a Latin Square, we try to fill the other slots of the CPLS with 
minimum number of symbols. 

From [7] it is known that if the Latin square L removes 
the singular fade state z then the Latin Square L T removes 
the singular fade state z _1 , where L T is the transpose of 
the Latin Square L, i.e. L T (i, j) = L(j,i) for all i,j E 
{0, 1,2, ..,M — 1}. This observation, in fact, holds for any 
choice of constellation S used at nodes A and B. 

From this, it is clear that we need to get Latin Squares only 
for those singular fade states with \z\ < 1 or \z\ > 1. 

The square QAM signal set has a symmetry which is tt/2 
degrees of rotation. This results in a reduction of the number of 
required Latin Squares by a factor 4 as shown in the following 
lemma. 

Lemma 4: If L is a Latin Square that removes a singular 
fade state z, then there exists a column permutation of L such 
that the permuted Latin Square L' removes the singular fade 
state ze J7T/2 . 

Proof: For the singular fade state z as given in |5]) with 
constraint {(xa, xb), (x' a , x' b )}, the singular fade state ze^l 2 



is given by 



oW2 _ 



W2 



[x' B - x B 
\x A ~ x' A ] 



A\ e J7r/2 



-J't/21 



Since in the square QAM constellation there exist signal points 
with x' B e~^/ 2 and x B e~° v ' 2 , though all the constraints are 
changed, the new constraints are obtainable by a permutation 
of signal points in the constellation used by node B. The 
columns of the Latin Squares are indexed by the signal points 
used by B and the effected permutation in the constellation is 
representable by column permutation in the Latin Square. ■ 
Lemma 5: If a Latin Square, L removes a singular fade 
state (7, 6) then the Latin Square to remove (7, (90 — 9)) is 
obtainable from L by appropriate row and column permuta- 
tions. 

Proof: See Appendix [B] ■ 
Note that the fade state z = 1 or (7 = 1, 6 = 0) is a singular 
fade state for any signal set. 

Definition 3: A Latin Square which removes the singular 
fade state z — 1 for a signal set is said to be a standard Latin 
Square for that signal set. 

It is known that the removal of the singular fade state z=\ 
has very high significance in a Rician fading scenario when the 
Rician factor K 7^ 0. The readers may refer to fl5) for further 
details regarding the influence of the values of singular fade 
states on the overall performance of the bi-directional relay 
network. 

When nodes A and B use a 2 A -PSK signal set, it has been 
shown in |7| that the Latin Square obtained by Exclusive-OR 
(XOR) is a standard Latin Square for any integer A. It turns 
out that for M-QAM signal sets the Latin Square given by 
bitwise Exclusive-OR (XOR) is not a standard Latin Square 
for any M > 4. This can be easily seen as follows: Any 
square il/-QAM signal set (M > 4) has points of the form 
a+jc, a+j(c+b), a+j(c—b), for some integers a, b and c. For 
2=1, the effective constellation at R during the MA phase 
contains the point 2(a + jc) and can be resulted in at least 
two different ways, since 2(a + jc) = (a + jc) + z(a + jc) = 
a+j(c+b) + z(a+j(c — b)) for z — 1. Let l\, I2 and 1% denote 
the labels for a+jc, a+j(c + b), and a + j(c—b) respectively. 
For the singular fade state z = 1, we have {(h, h), {h, h)} as 
a singularity removal constraint. But the Latin Square obtained 
by bitwise XOR mapping does not satisfy this constraint since 
h@h = 0^h®h- 
B. Standard Latin Square for y/M—PAM 

In this subsection, we obtain standard Latin Squares for 
a/M-PAM signal sets. 

Definition 4: An M x M Latin square in which each row 
is obtained by a left cyclic shift of the previous row is called 
a left-cyclic Latin Square. 

Lemma 6: For a \/M-PAM signal set a left-cyclic Latin 
Square removes the singular fade state z = 1. 



Proof: See Appendix [C] ■ 
Example 5: Consider the received constellation at the relay 
when the end nodes use 4-PAM constellation and let the 
channel condition be z = 1 as given in Fig [5] The singularity 
removal constraints are 

{(0,1)(1,0)}, {(0,2)(1,1)(2,0)}, {(0,3)(1,2)(2,1)(3,0)}, 
{(1,3)(2,2)(3,1)}, and{(2,3)(3,2)}. 

The Latin Square which removes this singular fade state is 
given in Fig|6] 

(0,1)(1,0) (0,3)(1,2)(2,1)(3,0) (2.3)(3,2) 
• • • • • • • 

(0,0) (0,2)(L1)(2,0) (1,3)(2,2)(3,1) (3,3) 

Fig. 5. Received Constellation at the relay for z = 1. 
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Fig. 6. Left-cyclic Latin Square to remove the singular fade state 2 = 1 

C. Standard Latin Square for M—QAM 

In this subsection standard Latin Square for a square M- 
QAM constellation is obtained from that of \AM-PAM con- 
stellation. 

Let PAM(i), for i = 1,2, ••• ,VM, denote the symbol 
set consisting of \[M symbols {(i — 1)VM, ((i — 1)\/M) + 
1, ((t - 1)%/M) + 2, •• • , ((i - 1)VM) + (VM - 1)}. Let 



LpAM{i) denote the standard Latin Square for VA/-PAM, 
with symbol set PAM(i) and also let Lqam denote the 
standard Latin Square for M-QAM. Then, Lqam is given 



in terms of L 



PAM(i)i 



1,2, 



>M, as the block left- 



cyclic Latin Square shown in Fig. [7] This is formally shown 
in the following Lemma. 

Lemma 7: Let PAM(i) for i = 1,2, • • • , vM, denote the 
symbol set consisting of \J~M symbols {(i — l)y/M,((i — 
l)v/M) + l,((i-l)>/M)+2,..- ,((i-1)Vm) + (VM-1)} 
and let L P am(i) stand for the Latin Square that removes the 
singular fade state z = 1 with the symbol set PAM(i) for 
y/M-PAM. Then arranging the cyclic Latin Squares LpAM{i) 
as shown Fig j7] where each row is a block wise left-cyclically 
shifted version of the previous row results in a Latin Square 
which removes the singular fade state z = 1 for M -QAM. 

Proof: See Appendix [P] ■ 
The standard Latin Square for 16-QAM is shown in Fig[8] 

IV. Channel quantization for M-QAM signal sets 



In Section III we saw how the relay R chooses a complex 
number for transmission during the broadcast phase if the 
channel fade coefficients result in a singular fade state. This 
was done essentially by associating a Latin Square to each 
singular fade state. However, 7e je being a ratio of two 
complex numbers (kg /Ha), can take any value in the complex 
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Fig. 8. Standard Latin Square Lqam f° r 16-QAM. 
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Fig. 7. Construction of L 



QAM 



for . 



plane (in fact, it takes a value equal to a singular fade state with 
probability zero). This demands the partition of the complex 
plane into regions and associating a Latin Square to each 
region so as to optimize total number of network coding maps 
used at the relay during the broadcast phase. The study of 
channel quantization for Af-PSK signal set for two-user fading 
MAC channel is described in p4| and that for bi-directional 
relay network is addressed in [8]. While [ 14] uses the notion of 
distance classes for partitioning the complex plane, in [ 8 1 the 
regions associated with each singular fade state are obtained 
by plotting the pair-wise transition boundaries of relevant pairs 
of singular fade states. 

In this section we describe how to partition the set of all 
possible ^e^ 9 when the end nods A and B use a square M- 
QAM signal set. The salient difference while using A/-QAM 
constellation when compared to A/-PSK signal sets at nodes 
A and B is as follows. It is shown in Lemma 1 of [9| that 
no two distinct pairs of points (d^, d/J and (d k2 ,di 2 ) in 
the difference constellation AS, such that \d kl \ ^ \dk 2 \ an d 



|eJ/J 7^ \di 2 \, can have the same ratio when S is Af-PSK signal 
set, i.e., 



d kl ^d h ,d k2 ^d h eAS (9) 



_ d kl _ dk 2 
d^ d h 
if and only if \d kl \ = \d k2 \ and \d h \ = \d h \. 



In other words, all z ^ z' ^ 1 € H, arise from distinct pairs 
differing in their magnitudes in AS for A/-PSK. This is not 
the case when the end nodes A and B use A/-QAM signal set. 
This is shown in the following example. First, notice that for 
A/-QAM signal set used at nodes A and B, any point d k € AS 
can be written in the form 



k = (n k + jm k ) , for some n k , m k G j- (VM - ij , 

-(-/M-2) J ...,(VM-2),(«/M-l)} 

(10) 



Example 6: Consider 16-QAM signal set used at nodes A 
and B. For the following 4 points in AS, d kl = -^=(2 + j) ^ 

dh =^(1 + 0j). d k 2 = 7=(1 + 3j) £ d l2 = ^L(l + j), 



I? 6 



d kl 
di t 



d^ 
di 3 



However, \d kl \ ^ |<4 2 | and \d h \ ^ \d h \. 

So the procedure used for channel quantization in (8) for 
A/-PSK does not directly apply for A/-QAM. However, as in 
case with A/-PSK described in |(8), we show that for certain 
values of ^e^ 8 any choice of clustering satisfying the exclusive 
law gives the same minimum cluster distance so that any one 
of the Latin Squares may be chosen. Subsequently, channel 
quantization for those values of je^ e for which the choice of 
the Latin Square does play a role in the overall performance 
is taken up. 



A. Clustering Independent Region 

Definition 5: The set of values of 7 e- je for which any clus- 
tering satisfying the exclusive law gives the same minimum 
cluster distance is referred to as the clustering independent re- 
gion. The region in the complex plane other than the clustering 
independent region is called the clustering dependent region. 

It is shown in [ 14] that with an arbitrary signal set S used at 
A and B for the MA phase, for any clustering C 1,9 satisfying 
the exclusive law, we have 



u inin 

(C 7 ' ) < mm{d m in(S), r yd mi n(S)} 



(11) 



Observation 1: Using ( fTTj ), for values of r ye J 
for which d min (je j6 ) > mm{d min (S), 7 d mm (S)}, 
since d min {^e je ) < d min (C^< e ), we have 

dminiC' 6 ) = mm{d min (S),id min (S)} for all clusterings 
C 1 ' 6 satisfying the exclusive law. Such je^ 6 therefore belong 
to the clustering independent region. 
Define 

T C i(S) ={~fe 36 : \d k +je je di \ > mm(d min (S), jd min (S)) 

V (4, d{) £ (0, 0) G (AS) 2 , 7 g R+ -it < 9 < ir} 

(12) 

(13) 
(14) 

ci ' ■ ' < '/ ' i - ' r • (15) 



where 



^}(S) u r%} (s) 



rg?(5) = r C z(5)n{ 7 >i} ; 



(<S) = r CJ (<S) n {7 < 1} • 



{4 



the set i± U/M - 
'M-lY- 



It has been shown in |8J that the region T™/ (<S) is obtained by 
the complex inversion of the region (<S) so that by finding 
out one the other can be obtained easily. 

Theorem 1: For A/-QAM signal set the region 
rgF/ (M-QAM) is the outer envelope region formed by 
— l] unit circles with centers (a + j/3) belonging to 

lj + jx, x ± j ys/M — ljX, where x G 

v / M-2),...,(VM-2),(n/M-i)} 

Proof: See Appendix \E\ ■ 
It can be verified that the centers of the 8 (^/M — 1^ circles 
in Theorem [T] are the singular fade states which lie on the 
outermost square. 

For M > 4 the region V™} (S) is described in the following 
lemma. Here, it may be noted that, for normalized signal sets 
used at nodes A and B, 4-PSK and 4-QAM are the same, and 
M-PSK (for all M) has already been addressed in (8). 

Lemma 8: For M > 4, rg^(M-QAM) is the 
outer envelope region formed by the 8(y/M — 1) 
circles with centers (- 

f3 2 — l) where (a + j/3) belongs 
M-lj +jx,x±j (yM-ljY x 

M-2V..., (Vm-2]JVm 



radii 1/ (a 2 4 
the set |± ( 

fJI - 1 



and 
to 

G 



{-6 



0} 



Proof: Proof follows directly from the complex in- 
version of the 8 {^\fM — lj circles used for obtaining 
rgF/ (M - QAM) in Theorem [l] ■ 



Define 



|| 7 e^ 



= max 



{M-ye?% |3( 7 e je )|} , V 7 e J ' e G C (16) 



where 3?(z) and ^s(z) stand for the real and imaginary parts 
respectively of the complex number z. 



From Theorem [l] it can be noted that for 



7 e je belongs to Tffj 



and for \ r )e> e \ < 



| 7 e- 



301 



> 



/237+1-V2' 



7 e- 



M, 
36 



je belongs to T % gj. In the complex plane, the locus of the points 



satisfying || 7 e 



3<? I 



a is a square centered at origin and 



having sides of length a. So, if 7 e je lies outside the square 
centered at origin and having sides of length \[M or inside the 

circle centered at the origin with radius r — 1 7= the relay 

can choose to use a fixed predetermined clustering satisfying 
the mutual exclusive law. In particular, that clustering whose 
corresponding Latin Square has only M symbols can be 
used. This observation helps in significantly reducing the 
computational complexity at the relay. The relay needs to 
adaptively switch between network coding maps only if the 
above two conditions are not satisfied. The following example 
illustrates Theorem Q] and the above observation. 

Example 7: Consider the case when the nodes A and B use 
16-QAM signal set. For this scenario, according to Theorem 
[l] the 24 unit circles centered at the singular states lying in the 
outermost square are sho wn in Fig. |9(a)| The region rg?/ is 



the shaded region in Fig. 9(a) The region T^f is the shaded 
region in Fig. |9(fe)| The circles in Fig. |9(6)| are those obtained 
by the complex inversion of the unit circles shown in Fig. 
9(a) The clustering independent region Tqi is the union of 



the shaded regions in Fig. 9(a) and Fig. |9(b)| Additionally, 
as an approximation, if 7 e Jt ' falls in the region between the 
square with side length 4 (Bl in Fig. [T0| and circle with radius 



3V2+1 



(B2 in Fig. 



maps. This is shown in Fig. 10 



10 1, the relay uses adaptive network coding 



B. Clustering Dependent Region 

In this section we partition that region of the complex plane 
where the choice of Latin Square significantly determines the 
performance of the bi-directional relay during the broadcast 
phase. It follows from Observation [l] that values of 7 e :,e for 
which d min (^e> B ) < mm{d min (S), jd min (S)} constitute the 
clustering dependent region. As the criterion described in 1 8 ] 
for partitioning the channel fade state complex plane when 
nodes A and B use A/-PSK depends on (|9]l which no longer 
holds for M-QAM as illustrated in Example d5l, we develop 
an alternative criterion for performing the channel partitioning 
in the clustering dependent region. 
With T> (7, 0, dk,di) defined as 

V{ 7) e,d k) di) = \d k +^ e d,\, (d k ,di) ^ (0,0) G (AS) 2 , 

(17) 

we have the following lemma. 




(a) Tgf/ (shaded region) for 16-QAM 




(b) Tg 1 / (shaded region) for 16-QAM 
Fig. 9. Cluster Independent regions for 16 QAM 



Lemma 9: If "fe J is such that 

org max ze ^ min (dfc di) . z £k =z 'D(j, 6, d k ,di) = z' , then 

the clustering maximizes the minimum cluster distance, 

among all the clusterings which belong to the set C-u 
Proof: For d k and di £ AS 

mm V{ 1 ,0,d kl d l )=d mm {C {z \ 1 ,0) 1 (18) 

(d k ,di):- s ^=z 

the minimum cluster distance of the clustering evaluated 
at je^ 6 . If the maximum in (18i is achieved for z' £ H, then 
d mm (d z '},j,e) > d mm {C {z \l,6) Vz ^ z' e H. Then 
is the best clustering for the channel fade coefficient 
je je . m 
Therefore, associated with each singular fade state z e T~L 
we have a region 1Z{ Z } in the je^ e plane in which the 
clustering C^ z ' maximizes the minimum cluster distance. This 
region 7?-{ z } is given by 



n 



{ 7 e 



i : min 

(d fc ,d,)G(A5) 2 



\d k + je je d t \ > 

— Z 

\d k +-fe je di\,yz' ytzeH}. 



The boundaries of these regions for each singular fade 
state are explicitly derived next. It is shown that like with 
A/-PSK signal set considered in [7], the boundaries of the 
region 1Z{ Z } are either circles or straight lines and a systematic 
procedure to obtain these regions for each singular state is 
given. A simple formulation to find out the pair wise transition 
boundary corresponding to a pair of clusterings C^ Zl ^ and 
C {z2} is stated next. 

The curve c(zi,Z2) which denotes the pair-wise transition 
boundary formed by the singular fade state z\ with the singular 
fade state z-i is the set of je^ 6 for which 

\d k +je ie d l \ = 



mm 

(d fc ,di)G(A5) 2 



mm 

(d k ,d,)G{AS) 2 



|d fc +7e ,w di|. 



Theorem 2: With the notations 




-5 -4 -3 -2 



-1 1 

Re(z) -> 



Fig. 10. The unshaded region shows the region in the complex plane where 
adaptive switching needs to be done for 16-QAM signal set 



di = arg min 
di = arg min 

d 2 £AS : -P~=2 



{|d 1 +7e J '°d a |} J 
{|di+7e J '%|}, 



the pair wise transition curve c (z\, Zq), z\ ^ Z2 is any one of 
the following 

• if \di\ ^ \df\, a circle with center (x, y) and radius r, 
where 



ft(*i) 
- 1^1 

' 1 di I 



ft(z 2 ) 



r, y 



l-l&l 



+ 



3(g) 



and r = 

if \d t \ = 
where 




Z2\ 2 \d,,\ 2 



2 Id/ 1 2 



|d/| 2 -|4 



|dj/|, a straight line of the form ax + by = c, 



a=(n(z 1 )\d l \ 2 -M(z 2 )\d l ,\ 2 ), 
b=(Z(z 1 )\d l \ 2 -Z(z 2 )\d l ,\ 2 ), 

c = ~(\z 2 \ 2 \d l ,\ 2 -\zi\ 2 \d l \ 2 ). 

Proof: See Appendix |f| ■ 
Observation 2: From (8| it is known that the region 
7Z{_di/d k } is the region obtained by the complex inversion 
of the region 7t{-d k /di}- Also the distribution of the sin- 
gular fade state in the ^e^ complex plane is periodic with 
periodicity tt/2 for M-QAM. Moreover, within an interval 
[a,a + tt/2) ,a <E {0, tt/2, it, 3tt/2}, the singular fade states 
are symmetric with respect to the line 9 = a + tt/A. 

Channel Quantization of 16-QAM Signal Set: We consider 
the case when nodes A and B both use 16-QAM signal set. 
From Observation [2] it follows that we need to consider only 
those singular fade states which lie outside the unit circle 
and within the angular interval 9 e [0,7r/4]. The pair wise 
transition boundaries as given by Theorem[2]for adjacent pairs 
of singular fade states in this region give the region TZ,{-d k /di}< 
\dk\ > \di\ and £{—dk/di) £ [0,tt/4]. Using the symmetry 
property and periodicity as mentioned in Observation [2] we 
get the region lZ{-d k /d,} corresponding to all singular fade 
states lying outside the unit circle centered at the origin. 
Again, from Observation [2] the regions corresponding to the 
singular fade states lying inside the unit circle is obtained by 
complex inversion of those obtained for singular fade states 
lying outside the unit circle. The regions corresponding to 
singular fade states lying on the unit circle are the remaining 
regions in the complex plane after the regions corresponding 
to all the other singular fade states are obtained. 

For 16-QAM signal set used at nodes A and B, there are 
27 singular fade states in the region outside the unit circle 
and in the angular interval 9 g [0, 7r/4]. To illustrate Theorem 
[2] consider the singular fade state 2 + j in the said interval. 
It shares pairwise transition boundaries with 8 neighbouring 
singular fade states viz. {1.5 + 1. 5j, 1.5 + 7, 1.4 + 0. 8j, 1.8 + 
0.6j, 2 + 2j, 3 + j, 2.5 + 0.5j, 2}. The pairwise transition 
boundary between 2 + j and 1.5 + 1. 5j is the circle CI shown 



in Fig. 11(a) Similarly, circles C2, C3, C4 and C5 form the 
pairwise boundaries with 1.5 + j, 1.4 + 0.8j, 1.8 + 0.6j and 
2.5+0.5j respectively, for the singular fade state 2+j as shown 
in Fig. 11(a) The pairwise transition boundary between 2 + j 



and 2 + 2j is the straight line LI shown in Fig. 11(a) Further, 



from Fig. 11(a) lines L2 and L3 are the pairwise transition 
boundaries with 3 + j and 2 respectively for the singular fade 
state 2 + j. Finally, the region TZ^+j} is the shaded region 



shown in Fig. 11(6) 



Proceeding likewise for all the singular fade states lying 
outside the unit circle, we get the channel quantization for 
16-QAM signal set as shown in Fig. [12] 

V. Simulation Results 

The Latin Square (LS) scheme (7) is based on removing 
the singular fade states. For 16-PSK all the 912 singular 
fade states can be removed with Latin Squares with number 
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(b) Region K{ 2 +j} for 16-QAM 

Fig. 11. Diagram explaining the procedure to get region for a singular fade 
state 

of symbols 16, but for 16-QAM some singular fade states 
cannot be removed with Latin Squares of 16 symbols, we 
used Latin Squares consisting 20 symbols for some singular 
fade states. Since 16-QAM has only 388 singular fade states, in 
comparison with 912 singular fade states of 16-PSK, and since 
16-QAM offers better distance distribution in the MA stage 16- 
QAM gives better performance. For a given average energy, 
the end to end BER is a function of distance distribution of 
the constellations used at the end nodes as well as at the relay. 
The simulation results for the end to end BER as a function 
of SNR is presented for different fading scenarios. 

Consider the case when hA,h,B,h' A and h! B are distributed 
according to Rayleigh distribution, with the variances of all the 
fading links are assumed to be dB. The end to end BER as a 
function of SNR in dB when the end nodes use 16-QAM signal 
sets as well as 16-PSK signal sets with same average energy 
is given in Fig[l3] The end to end BER for XOR network 
code for 16-QAM is also given. It can be observed that the 
LS Scheme for 16-QAM outperforms LS Scheme for 16-PSK 
as well as XOR network code. 




■MS Scheme for 16-QAH 
-«-LS Scheme for 16-PSK 
-•-XOR NW Code for 16-QAM 
XOE NW Code for 16-PSK 



Fig. 12. Diagram showing the regions associated with the singular fade states 
Wins outside the unit circle for 16-OAM sienal set 



-•-LS Scheme for 16-QAM 
♦II Scheme for 16-PSK 
-»-X0R NW Code for 16-QAM 
XOR NW Code for 16-PSK 




Fig. 13. SNR vs BER for different schemes when the end nodes use 16-QAM 
and 16-PSK for a Rayleigh fading scenario. 



Consider the case when hA,h,B, h' A and h' B are distributed 
according to Rician distribution, with the Rician factor of 5 
dB and the variances of all the fading links are assumed to be 
dB. In Fig{l4]the end to end BER as a function of SNR in 
dB for LS scheme for 16-PSK, 16-QAM and XOR network 
coding for 16-QAM is given. It is observed that the LS scheme 
gives large gain over the XOR network coding scheme. The 
LS scheme for QAM is better in end to end BER performance 
in comparison with the LS scheme for PSK. 

VI. Conclusion 

In this work, the design of modulation schemes for the 
physical layer network-coded two way relaying scenario when 
the end nodes use square QAM constellation is studied. We 
show that there are many advantages of using square QAM 
constellation over PSK signal set. Construction of the standard 
Latin square for removing the singular fade state z=\ for M- 
QAM is described and this is shown to be different from that 
used for M-PSK. Using the relation between exclusive law 




Fig. 14. SNR vs BER for different schemes when the end nodes use 16-QAM 
and 16-PSK for a Rician fading scenario with Rician factor 5 dB. 



satisfying clusterings and Latin Squares, a method to remove 
all the other singular fade states is proposed and a means to 
derive the corresponding Latin squares is presented. This gives 
us all the maps to be used at the relay when square QAM 
constellation is used at the end nodes. The channel partition 
for QAM signal set is obtained analytically. Simulation results 
showing the end to end BER performance when the end nodes 
use PSK constellation as well as QAM constellations are 
obtained to support our claim. 
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Appendix A 
PROOF OF LEMMAQ] 

There are 2{yM — 1) non-zero signal points in the differ- 
ence constellation AS and since AS is symmetric about zero 
there are yM— 1 signal points in AS + . All these are scaled 
version of nonzero elements of Z v /jy. 

The number of positive integers less than or equal to n that 
are relatively prime to n is given by, 

^(n) = njj(l-i) (19) 

p\n 



P 



where the product is taken over distinct prime numbers p 
dividing n. To get the total number of relatively prime pairs in 



M 



, we take the sum over all nonzero neZ fjj which gives 



1 ^n p |n( 1 -^)- ° ne relatively prime pair (a, b) 
gives two singular fade states, a/b and b/a. The multiplication 
factor 4 in |6]) accounts for the negative side of the in-phase 
axis as well as the inverses. Finally, 2 is added to count the 
singular fade state z = 1 and z = — 1. ' 



Let 7e J 



tan 



-1(JQ_\ _ 

V 7I / 



Appendix B 
PROOF OF LEMMAS 

7/ + J7Q» such that 7? + 
9. Then 7 e^ 9 °- e ) = 7q + jjr. We have, 



and 



[XA 



X, 



Let xa ■ 
Then, 



Xd 



7e j« = 

[X'b - X B \ 

j+jxd A Q and x' B - 



xb = Xd B I+jXd B Q. 



7/ + 3lQ 



Xd, 
Xd, 



JXd A Q 
jXd B Q 



After expansion, we get 

{Xd A lXd E 



II + 31Q 



I + X dA QXd B Q, 



+ 



and 

1Q 



■311 



j{x dA QXd B I - X dA lXd B Q) 

( X d B I + x d B o) 



{Xd A QX dB I - X dA lX dB Q) 

( X d B I + x d B Q) 

j(x dA lX dB I + X dA QX dBQ ) 

T „,2 n 



(20) 



(21) 



By comparing above expressions we can say that Xd A i, 
Xd A Q ,Xd B i and x dB Q in (|20]i are changed to x dA Q, x dA i, 
Xd B i and —Xd B Q in ( pi) , i.e., the difference constellation 
points x dA i + jx dA Q and x dB i + j x d B Q whose ratio gives 
singular fade state (7, 9) are converted to Xd A q + jxa A i 
and Xd B i — jxd B Q and whose ratio gives singular fade 
state (7, (90 — 9)). Let Xk — x^i + jxkQ, where k e 
{A,B,A',B'} and let the singularity removal constraint 
for ( 7 ,0),be {xa,x b ){x' a ,x' b ) = (x A i + jx A Q,x B i + 
JXbq)(x' ai + jx' AQ ,x' BI + jx' BQ ). The singularity removal 
constraint for (7, (90 — 9)), can be written as (x A q + 
jx A i,x B i - JXbq)(x' aq + jx' AI , x' BI - jx' BQ ). Since con- 
stellation points used by node A are indexed by rows in the 
Latin Square, the rows corresponding to constellation point 
x A i + jx A Q are permuted to the rows denoted by xaq +jxAi 
and since the signal points used by node B are indexed by 
columns in the Latin Square, the columns corresponding to 
constellation point xbi + ]Xbq are permuted to the columns 
denoted by xbi — JXbq, to get the Latin Square to remove 
(7,(90-0)). ' ' " ■ 

Appendix C 
PROOF OF LEMMA [6] 

Consider the ^/M-PAM signal set with the signal points 
labelled from left to right as discussed in Section [TTJ Let 
{(&!, h)(k2, h)} be a singularity removal constraint. To get 
the same point in the received constellation at the relay R, 
when z = 1, we have fci+Zi = k2 + h- Consider the following 
two cases satisfying this equality. Case (i): /C2 = h,h = ki 
In this case the constraint becomes {(/ci, Zi)(Zi, fci)}, i.e., the 
Latin Square which removes z = 1 should be symmetric about 
main diagonal. 

Case (ii): &2 = &i + m, h = h ~ m for any m < \/M., 
The constraint now becomes {(fci, l\){ki +m, l\ — to)} which 
means the symbol in fci-th row and Zi-th column should be 
repeated in the k\ + 1-th row and the l\ — 1-th column. 
It is easily seen that a left-cyclic Latin Square satisfies both 
this requirements. ■ 

Appendix D 
PROOF OF LEMMA [7] 

Note that the matrix in Fig. [7] is a M x M matrix, which is 
also a \[M x \J M block left-cyclic matrix where each block 



P AM (i) 



for some i. 



is a y/M x \JM left-cyclic matrix L 

Let ai+jbi,a2+jb 2 ,a' 1 +jb' 1 and a 2 +jb 2 , where aj,a^,6j 



and b[ e {-(VM-l),-(VM-3),— , (VM - 3), (VM - 
1)} for i E {1,2} be four M-QAM constellation points such 
that <zi + jbi and a[ + jb[ are used by node A and a 2 + jb 2 
and a' 2 + jb' 2 are used by end node B, and result in a same 
point in the effective received constellation at the relay node 
for singular fade state z = 1, i.e., 



ai + jh +a 2 + jb 2 = a\ + jb' t 



■3b' 



Let a[ - 

{-2(VM- 



f toi and b[ 

-2(VM-2),- 



h - 
,2( 



to 2 where toi , TO2 e 
M-2),2(\/M-l)}. 



Then, a' 2 = a 2 — mi and b' 2 = b 2 — m 2 . Then, using the map 
defined in {7]), let 

h = t4 a i + 3 b i), 
h = m(«2 +jh), 
k 2 = li{a\ + jb[) = /i(ai + mi + j{b\ + m 2 )) and 
h = li{a! 2 + jb 2 ) = fi(a 2 -mi+ j(b 2 - m 2 )). 

Since, for z = 1, the four complex numbers result in 
the same point in the effective constellation at the relay, 
{(fci, l\)(k2, h)} is a singularity removal constraint for z = 1. 
From the above equations it follows that 

1 



k 2 = ki + - (miVM + m 2 ) 



h = h- 2^ mr 



^ + m 2 ), 



The above equations precisely mean the construction shown 
in Fig|7] This completes the proof. ■ 

Appendix E 
PROOF OF THEOREMQ] 

From (jXSJ and ( fl4] i we have 

Ic?(5) = {ie ]6 ■ \d k + l^ 9 di\ > min(d mi „(<S), 7 d min (5)) 

V(rf fc ,d;) 7^ (0,0) € A«S X A«S,7> 1,-tt< 6< tt}. 

(22) 



Let ~d k /di 
we have 



z€?i with |z| > 1. Then, Vje 3 ' 6 G T^f/ (5), 



K| 2 |-z + 7 e J T > (d min (S)f 
Vd ; ^ G A5,Vz = - 



-dfe/di G H with |z| > 1. 

(23) 



In particular, M^e ie G rgf (5) 

{|d ; | 2 |-^ + 7e J T} > (4 



mm 



(«S))' 



(24) 



Now min dl ^06A5 {Mi| 2 } = (d mm (S)) 2 = 4/p, from 

ffor M-QAM. Using <[ToJ, \d t \ = implies that 

is of the form ±2/yfp or ±j2/y/p. For such d;, 
z = —dk/di is of the form a + jf3, where a,/3 G 5" = 
M - lV - f VM - 2^ , (VM - 2 s ) , L/M - l" 



(25) 



so (|24li reduces to the form 

|- a- j/3 + 7e 3 T > 1, Va,^e5". 



This is the exterior of the unit circle with center 
(a, P). It must be noted here that VjeJ 8 G 
rgy*(5) and Va,/3 



(25 1 holds for the outer envelop region of the unit circles 
with centers (a, f3). This means 

a + i^e{±(VM-l) +jx,x±j(VM-l)}, 

where (26) 



From (26 1 the number of such circles is 2 x 2 ^2%/M — 1 



'M - 1 . 



Appendix F 
PROOF OF THEOREM[2] 

At the pair wise transition boundary corresponding to the 
pair of clusterings c{~ dfc i/ dl i} and c{ _d *2/*2 J, denoting 
—dki/d^ as z\ and ~d k2 /di 2 as z 2 , the set of values of je 3 
satisfies 

min {\d kl +je je di 1 \} = min {\d k2 + ^e> 6 di 2 1} . 

d^: di 2 : 

(27) 

Say the min in the LHS and RHS of (27i are achieved at 



di and d;/ respectively. Using the notation used in (lOi, (27 1 
becomes, 

l^fei + jm kl + je 30 (n t + jrhi) | = 

\n k2 + jm k2 + je 3e {ri v + jra v ) |. (28) 



With je 3 — 7^ + j-fi, squaring and rearranging (28i gives 



I - 2 i ~ 2 ~ 2 - 2\ 2 , / - 2 , -2 -2 - 2\ 2 



+ 2 (n fcl n ( + m kl mi - n k2 ri v - m k2 m v )^ R 
+ 2 {m, kl ni - n kl rhi - m k2 ri v + n k2 vH v )^i 



Now, since 



n k2 + m k2 - n kl 2 - m kl 2 



n fe n ; + m k mi 



(29) 



(30) 



and 



9 (-*/*) = 9 (-=4^) 
V ni+jmij 

_ m k ni - n k mi 
nf + m 2 



(31) 



29j) becomes 

(|d ; | 2 - | 2 ) ( 7 i? 2 +7/ 2 ) - 2 (R(*i) |d ; | 2 M"/'| 2 ) 
- 2 |d ; | 2 - 3 (z 2 ) \d v \ 2 ) = \z 2 \ 2 \d v \ 2 - M 2 K| 2 

(32) 



G 5" (l25l) must hold If ^ M'l 7^ Mi' I > < 32 » is the equation of a circle of the form 



{a,0)€ [S»-{±{VM-1)}) x 5"-{±(VM-l)} , 



<|25b does not hold for (a,/3) £ {±(VM-l\\ x 
■ ± (^/M — So either one or both of a and /3 must 
belong to the set |± (^/M — l\\ and for such a choice, 



(th ^ a: ) 2 + (7/ - y) 2 = 7,2 

with x, y and r as given in Theorem [2] 

If \di\ 2 = |d/'| 2 , (32 1 is a linear equation. Since z\ 7^ z 2 , the 



coefficients of both 7^ and 7/ cannot be simultaneously zero, 
so that we get the equation of a straight line. ' 



